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We determine the classical and quantum complexities of a specific ensemble of three-satisfiability problems 
with a unique satisfying assignment for up to N = 100 and TV = 80 variables, respectively. In the classical 
limit we employ generalized ensemble techniques and measure the time that a Markovian Monte Carlo process 
spends in searching classical ground states. In the quantum limit we determine the maximum finite correlation 
length along a quantum adiabatic trajectory determined by the linear sweep of the adiabatic control parameter 
in the Hamiltonian composed of the problem Hamiltonian and the constant transverse field Hamiltonian. In the 
median of our ensemble both complexities diverge exponentially with the number of variables. Hence, stan- 
dard, conventional adiabatic quantum computation fails to reduce the computational complexity to polynomial. 
Moreover, the growth-rate constant in the quantum limit is 3.8 times as large as the one in the classical limit, 
making classical fluctuations more beneficial than quantum fluctuations in ground-sate searches. 

PACS numbers: 75.10.Nr, 02.70.Ss, 03.67.Ac, 64.70.Tg 
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At the borderline of mathematics and physics lie optimiza- 
tion problems, that can be cast into solving a minimization 
problem on a discrete set of variables: Given a scalar cost 
function Hq(s) that is bounded from below by zero, and 
given a set of integer variables (Ising spins) s, = ±1 with 
i = 1,...,N, one may ask: Which assignment or - satisfy- 
ing assignment - solves Hq =01 Many satisfying assign- 
ment problems with Boolean variables b{ = (1 + Sj)/2 that 
are NP-hard have precisely this form. In this work we study 
the 3-satisfiability (3-SAT) problem, a NP-hard problem at the 
heart of complexity theory yl, by means of methods used in 
(quantum) statistical physics. 

Under the assumption P ^ NP, the computational effort 
for any classical algorithm to solve NP-hard problems is be- 
lieved to be 0(e gN ), where g denotes the growth-rate con- 
stant. For a trivial classical or unstructured search, that is the 
evaluation of Hq over all configurations, g — In 2. In quantum 
computation this search finds its analog in Graver's algorithm 
Q2|] with g = In 2/2. A polynomial solution to a NP-hard 
problem is expected to have g = 0. 

Conventional, standard adiabatic quantum computation 
(AQC) JH 0] assumes a linear interpolation between the NP- 
hard problem Hamiltonian Hq and a non-commuting driver 
Hamiltonian Ho = J^i a i ( me "transverse field"), where erf 
is the x-component of the Pauli matrix. A statistical analysis 
of AQCs determines the thermodynamic and quantum singu- 
larities of the partition function 

Z AQC (l3,\)=Tre-W l -» H ° +XH °\ 



(1) 



where /3 = 1 jk^T denotes the inverse temperature, ks Boltz- 
mann's constant and < A < 1 the quantum adiabatic control 
parameter. In the vicinity of the point Pq* = (f3, A) = (oo, 1) 
the optimization problem is solved, as vanishing thermal and 
quantum fluctuations lead to the exact ground-state. We study 
the approach to Pq from regions of large thermal, as well as 



large quantum fluctuations on lines of parameters (3 and A. In 
particular, we study measures of complexity in the classical 
limit at A = 1 as a function of (3 and in the pure quantum limit 
at almost zero temperature as a function of A. 

In the classical limit, a measure of complexity is the Monte 
Carlo (MC) search time for the ground state in multicanonical 
ensemble and Wang-Landau simulations J5[0]. These MC 
simulations perform a Markovian process with random walk 
dynamics in the energy. We count the number of MC steps 
in-between ground-state findings in the mean. 

In the pure quantum limit we determine the maximal spin- 
spin correlation length £ max , i.e. the inverse of the first energy 
gap at the presumed quantum phase transition at A*, from the 
exponential decay of a two point function in imaginary time 
(see Eq. ©. If there exists an avoided level crossing, the spin- 
spin correlation length is finite for a finite number of spins. In 
accordance with Landau Zener theory H, in AQCs the run- 
ning time of ground-state searches is limited to a time scale T 
of order 0(£max) fr° m below. Hence, for a NP-hard problem, 
a spin-spin correlation length growing exponentially with N 
would yield a computational complexity for quantum ground- 
state searches that is similar to the one expected for a classical 
search, and therefore would make AQC fail. 

It is argued that exponentially small energy gaps can be in- 
duced by the presence of first-order phase transitions, ham- 
pering the performance of AQC in optimization problems re- 
lated to the 3-SAT problem Il8l 4llll . Their predictive power 
for specific optimization problems is however limited, as any 
first order phase transition may - or just may not - turn into 
second order at a critical point. A particular nasty situation is 
encountered, if there is either a weak first - or second order 
phase transition, a situation that has recently been studied in 
the exact cover problem [12]. 

The 3-SAT problem is defined on a set of i — 1, 7Y clas- 
sical Ising spins Sj = ±1. Its Hamiltonian can be written as a 
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FIG. 1 : Classical limit: Expectation value of the problem Hamilto- 
nian {Ho) as a function of /3 for a particular realization with N — 80 
and r — 8. Arrow: Position of the maximal specific heat. Inset: Dis- 
tribution function of the ground-state overlap observable, see Eq. 10). 
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FIG. 2: Classical limit: Expectation values for the logarithmized 
search time \n{T ( J /N 2 ) (triangles) and the nucleation barrier {Bo) 
(circles) as a function of the number of spins N. Straight lines corre- 
spond to exponential singularities. Open symbols: r — 8; full sym- 
bols: r = 5. Inset: N = 100 and r = 8 dependency of ln[r°/N 2 } 
on Bo within a set of 1000 realizations. 



sum of M three-point functions, called clauses: 
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Tf 3 (e aj iS i [ Qi i], e a ,2Si[ Q ,2], e<x,3Si[a,3])- (2) 



The function T 3 results from a transcription of the disjunctive 
cubic Boolean form of a clause bi V b k V b m to Ising degrees 
of freedom: 

T 3 (Sfc, Sl, S m ) = -{(s k SlS m ) + (SfeS/ + SiS m + S k S m ) 
o 

+{s k + .n + s m ) - 1}, (3) 

with k,l,m e {1,...,N} and k ^ I ± m. If p = 1, 2, 3 
denotes the position of a spin within clause a and if i] labels 
an element from a set of TV^ realizations of spins Sj, the quan- 
tity p takes values ±1. The symbol i[a,p] denotes a map 
[a, p] —± i from indices of the clause to the set of spins. 

The 3-SAT problem has been studied on random instances 
with statistical methods, the hardness being characterized by 
the clauses-to-variables ratio r — M/N. It exhibits at least 
two phases: a "SAT" phase for r < 4.2 with degenerate 
ground-states and an "UN-SAT" phase for r > 4.2 where sat- 
isfying assignments are exponentially rare lf]~3h . Note that also 
computationally demanding instances can be constructed for 
r = 3 H. 

We study instances with unique satisfying assignments 
(USA) and for which r is largely free. For this pur- 
pose we generate random ground-states sf s , and solve 
£a,isf[ S Q i] = 1 for a particular map i [a, 1] = mod(a— 1, N) + 
1 and for a = 1,...,M. The remaining e a , p and i[a,f>] for 
p = 2, 3 are generated at random with the help of MC up- 
dates. We use various filter techniques to exclude non-USA 
instances. 



In studies of the classical limit for r = 5, 8 we use mul- 
ticanonical ensemble and Wang-Landau simulations 0] 
and determine a statistical estimate of the density of states 
n(E) on the entire discrete energy interval < E < E max . 
We obtain the canonical partition function Z can (j3) via its 
spectral representation and calculate thermodynamic quanti- 
ties like the internal energy 9 ( 3lnZ can (/3) and the specific heat 
/3 _2 i9|lnZ can (/3). Numerical analysis shows a discontinuous 
phase transition at some value /3*, see Fig.Q] Its first order 
nature can best be established by considering the ground-state 
overlap observable o gs = A^ 1 s i s T an d i ts distribution 

function at the specific heat peak position /3cmax ~ P* 



P( 0gs ) = Z- 1 ^'<5 



1 N 

i=l 



e" &MAxB , (4) 



where Y^, denotes a sum over all spin configurations. For 
almost all realizations P(o gs ) exhibits a bimodal shape with 
one sharp peak at o gs , r i g ht ~ 1 with a value P m ax, right; an ex- 
ample being depicted in the inset of Fig. Q] The second peak 
with value P raax ,ieft i s we U separated from the first one and 
is located at o gS) i e ft. At r = 8 a finite size scaling analysis 
yields a non-vanishing overlap order parameter gap (Ao gs ) = 
(ogs.right) — (ogs.icft) = l-H(l) for the thermodynamic limit 
and in the mean of realizations. Discontinuous phase transi- 
tions are also associated with nucleation free energy barriers. 
As the maxima in P(o gs ) are separated by a minimum with 
value P m i n , and in analogy to free energy barrier definitions 
in Ising magnets and glasses, we use Binder's method [ 15] to 
define a nucleation barrier B = ln[P max i c f t /P min ] 
for each realization r). 

We also perform MC simulations in the multicanonical par 
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tition function Z muca = ^ exp[— In n(E)] yielding an equal 
distribution for the probability P m uca(^) = const to find an 
energy E in the Markov chain. The MC dynamics in the 
energy E is different from that of a free random walk with 
polynomial singularities r oc iV 2 in autocorrelation times, as 
hidden free energy barriers at energies in the vicinity of the 
ground-state E = slow down the diffusion. Realization 
wise we measure the number of local Metropolis update steps 
Tg /N 2 in the mean, that an interacting walk in energy spends 
in the "transition" from E max towards the ground-state energy. 
The factor N~ 2 corrects for the trivial free walk behavior that 
is present even in the absence of barriers. 

Figure [2] displays the results for the averages of ln[r°/iV 2 ] 
(triangles) and Bq (circles) calculated in the median for r — 
5, 8 and a number of spins up to N — 100. Both quantities 
exhibit exponential behavior of the form A exp[g c N] with 
classical growth-rate constants g r c = 0.077(4) (barrier) and 
= 0.078(3) (search time) for r = 5. For r = 8 we find 
g T c = 0.016(1) and gf = 0.021(2), respectively . Hence, as 
expected, the complexity grows as r is lowered to r* s» 4.2 
from above. Note however that the growth-rates are much 
smaller than g — In 2 of an unstructured search. The inset of 
Fig. [2] displays binned mean values of r° /N 2 as a function 
of the nucleation barrier Bq for r = 8 and N = 100 for a 
set of N v = 1000 realizations. Within this set of realizations 
violent fluctuations of complexity related observables are ob- 
served. The search time also shows an exponential behavior 
of the form t°/N 2 = a + b exp[oB ] with c « 0.8. Thus the 
low-temperature free energy landscape of the 3-SAT problem 
has the simple property that the static free energy barrier Bo 
determines the ground-state search dynamics. 

We now quantize the problem by introducing a standard 
Trotter-Suzuki time discretization II 1611 . We choose a regular 
temporal lattice with iV r = 128, 256 time-slices, a finite step- 
size At = 1 in r direction and periodic boundary conditions 
in Trotter time. The inverse temperature is j3 = N t At and 
the Boltzmann factor of the quantized problem at imaginary 
time is 

N T 

ln[P B , q ] = -Ko^Hods^r), s n (t)}) 



N N T 



(5) 



with positive ferromagnetic hopping parameters kq — AAt 
and k t = — ln[tanh((l — A)Ar)]/2 Biol . These equations 
implement the AQC partition function Eq. (HJ as a function of 
A and /?, up to discretization errors caused by the finiteness of 
the regularization. The first energy gap is obtained from 
a large r fit to the exponential decay of the expectation value 
of the connected two point function in the canonical mean 



r(r) = (O(O)O(r)) - (O(0)) 2 rx e 



(6) 




We have experimented with several time-local observables 
0(t) and found that among an extended set of trial ob- 
servables the r-dependent ground-state overlap Cgs(r) = 



FIG. 3: Quantum limit: Signature of the quantum phase transition for 
the same realization of Fig.fTJ The mean overlap observable (o gs ) as 
a function of A shows a jump at A* ~ 0.446. Inset: Overlap-overlap 
correlation function T(r) at A* as a function of imaginary time r. 



TV -1 Yli=i s i( T ) s T yields the best statistical signals for the 
exponential decay with r. The minimum energy gap, or the 
maximum spin-spin correlation length £ max determines the 
complexity. We find it by a search, that uses parallel temper- 
ing simulations in the control parameter A IU7I1 on a A partition 
with 64 elements. Elementary low-temperature exchange up- 
dates are absolutely essential for error reduction in the quan- 
tum correlator. 

We have studied the quantum phase transition of the parti- 
tion function Eq. ([1) for r = 8 and N values up to N = 80. 
In Fig. [3] we display the expectation value of the quantum 
ground-state overlap observable (o gs ) = N~ 1 J2n °gs( T ) 
for a specific N = 80 realization. A quantum phase transition, 
which is of blatant discontinuous nature, is observed at A*. 
Interestingly, the transition to the ground-state proceeds with 
increasing A from configurations that are slightly anti-parallel 
to the known ground state, see the negative values of o gs for 
A < A*. The inset of Fig.[3]displays the quantum correlation 
function T(r) of Eq. © at A*, which for the specific example 
decays exponentially at the numerical value £ max « 185. 

In the analysis of the quantum complexity we have picked 
a subset of 91 realizations for N = 60 at about median - 
or less - classical complexity given by the Bq value. In the 
inset of Fig. [4] we display the correlation of Bq values with 
their quantum counterparts, the maximal spin-spin correlation 
length £ max values. We observe a linear correlation between 
m £max and Bq, see the straight line, and not even the slightest 
tendency that the quantum complexity is weaker than the cor- 
responding classical one. A similar observation is made for 
all N. 

We simplify the calculation of the quantum complexity by 
use of the median. Note that median averages do not necessar- 
ily require the actual calculation of extremal values, but just 
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FIG. 4: Quantum limit: Finite-size scaling of ln(£) as a function of 
N for r — 8. The largest median correlation length is (£) = 259.9 at 
N = 80. The straight line displays the fit Eq. Q7). Inset: Scatter plot 
of tupels of In £ and Bq (triangles) for N = 60. The data exhibit a 
linear correlation (straight line). We determine the correlation length 
at the position of the median barrier (Bo) = 1.90(6) (arrow and 
solid circle). 



of those in the vicinity of the median. Given our data, we may 
determine the spin-spin correlation length £ at the median po- 
sition of the nucleation barrier (-Bo)- For correlated data the 
latter quantity is a good estimator of the spin-spin correlation 
length in the median (£). A fit to the form 



(0 = A 



N 



(7) 



on the entire iV-interval 20 < N < 80 with Xd.o.f = 0.1 
yields for the quantum growth rate constant g q — 0.061(1) 
for the 3-SAT problem with r = 8. Figure|4]shows ln(£) as a 
function of N (circles) and the fit in Eq. (Q (straight line). As 
can be clearly seen, the data are not compatible with a polyno- 
mial behavior. Following the Landau-Zener theory, twice the 
growth rate constant 2g q should be considered for the compar- 
ison to the classical search time complexity. 

In summary, we have determined the exponential singulari- 
ties that dominate the classical and quantum running times for 
ground-state searches in the median of an ensemble in the 3- 
SAT problem with unique satisfiability assignments. For r = 
8 we measure the growth rate constants g c = 0.016(1) and 
g q = 0.061(2) of the corresponding classical t° and quan- 
tum (, ma x exponential behavior in the ground-state searches. 
Our numerical data are precise and the classical growth rate 
constant is confirmed by a static free energy i.e., nucleation 
barrier Bo scaling. The finite size scaling window is large 
and excludes a polynomial behavior for quantum and clas- 
sical searches. For the case of standard AQC we find that 
twice the growth rate constant 2g q undershoots the Grover 
value g — In 2/2. Hence, standard AQCs constitute a class 
of ground-state searches, that compares favorable to Grovers's 



quantum search. However, there is no indication that for the 
3-SAT problem standard AQCs can reduce the exponential 
complexity to a polynomial one. Thus, standard AQC ground- 
state searches for NP-hard problems are not quite as powerful 
111 811 as was conjectured earlier [4]. We also find a quantum 
growth rate constant that exceeds the classical one by a factor 
of 3.8. Therefore, sole quantum fluctuations are even signif- 
icantly less efficient than sole classical fluctuations. In the 
future one may determine complexities of systems with mod- 
ified driver Hamiltonians, that are optimized with respect to 
ground state searches. An alternative quantitative approach to 
computational intractability can rely on free energy calcula- 
tions for Instantons [19]. 
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